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Abstract 

The statefinder parameters (r, s) in two dark energy models are 
studied. In the first, we discuss in four-dimensional General Relativity 
a two fluid model, in which dark energy and dark matter are allowed 
to interact with each other. In the second model, we consider the DGP 
brane model generalized by taking a possible energy exchange between 
the brane and the bulk into account. We determine the values of the 
statefinder parameters that correspond to the unique attractor of the 
system at hand. Furthermore, we produce plots in which we show s, r 
as functions of red-shift, and the (s — r) plane for each model. 



1 Introduction 

A plethora of observational data are now available, which show that we live in 
a flat universe that expands with an accelerating rate and that the dominant 
component in the energy budget of the universe is an unusual material, the 
nature of which still remains unknown. Identifying the origin and nature of 
dark energy is one of the great challenges in modern theoretical cosmology. 
The simplest candidate for dark energy is the cosmo logical constant, which 
corresponds to a perfect fluid with state parameter w = p/ p = —1. The 
LCDM model is still in agreement with all observational data. However, due 
to the problems associated to the cosmological constant, over the years many 
other theoretical models have been proposed and studied. One class of such 
models is based on some modification of Einstein's gravity [H [2] and one is 
talking about the so-called geometrical dark energy models. Another class 
contains the dynamical dark energy models, in which a new dynamical field 
(almost certainly a scalar field) is coupled to gravity. In the second class 
one would find models called quintessence [3], phantom [I], quintom [5], k- 
essence [Bj, tachyonic [7] etc. A recent review on dark energy dynamics one 
can find in [Sj. Furthermore, an attempt to solve, or at least to alleviate, the 
coincidence problem is realized in models of interacting dark energy [9j. In 
this class of models an interaction between the dark matter and dark energy 
is allowed, which after all is a natural possibility to be considered. 

Nowadays several cosmological models of dark energy are available which 
cannot be excluded by current observational data. However, a few years 
ago two new cosmological parameters were introduced [10] in order to dis- 
criminate between different dark energy models. These are the so called 
statefinder parameters and are given entirely in terms of the scale factor and 
its derivatives with respect to the cosmic time, up to the third order. The 
statefinder parameters are defined as follows 



aH 3 
r — 



(2) 



where a is the scale factor of the universe, H is the Hubble parameter, a 
dot denotes differentiation with respect to the cosmic time t and q is the 
deceleration parameter 



After the deceleration parameter with the second derivative of the scale factor 
with respect to the cosmic time, a natural quantity to be considered is r with 
the third derivative of a with respect to t. For the LCDM model, r = 1 and 
s = 0. At this point the easiest way to see this is to consider the evolution 
of the universe at large times, that is when the contribution from matter is 
negligible. Then the universe expands due to the cosmological constant only 
and the scale factor grows exponentially with time, a(t) = exp(t). Then one 
can immediately see that r = 1 and therefore s = (since q=-l). 

The trajectories in the r — s plane for various existing models can ex- 
hibit quite different behaviors. The deviation of these trajectories from the 
(0, 1) point defines the distance of a given model from the LCDM model. 
The statefinder pair (r, s) can successfully differentiate between a wide vari- 
ety of cosmological models including a cosmological constant, brane models, 
quintessence, Chaplygin gas, and interacting dark energy models. In a given 
model the pair r, s can be computed and the trajectory in the r — s plane 
can be drown. Furthermore, the values of r, s can be extracted from future 
observations [TT]. Therefore, the statefinder diagnostic combined with future 
observations may possibly be used to discriminate between different dark 
energy models. 

Up to now, the statefinder diagnostic has been applied to several models, 
see e.g. [12J. In the present work we wish to study two dark energy mod- 
els. In the first we consider in four-dimensional General Relativity a two 
fluid model with dark matter and dark energy interacting with each other. 
The form of the interaction is specified below. We do not rely on a con- 
crete particle physics model for dark energy. We just treat dark energy as 
a hydrodynamical fluid with a constant state parameter w = px/ Px, where 
Px,Px are the pressure and energy density of dark energy respectively. In 
the second model, we consider the DGP brane model generalized by taking 
into account a possible energy exchange between the brane and the bulk. 
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2 The dark energy models 



2.1 Interacting model 

The equations of motion for our system are Friedmann equations and the 
semi-conservation equation for each fluid component 

H 2 = yP (4) 

H = -y(p + p) (5) 
Q = p m + 3Hp m (6) 
-Q = Px + 3H(p x + Px ) (7) 

where Q is a source term responsible for the interaction between the two 
fluid components, k 2 = 8nG, p m is the energy density of matter (p m = 0), 
P = Pm + Px is the total energy density and p = p m + px = Px is the total 
pressure. The first Friedmann equation is a constraint, while the second one 
is a dynamical equation. Below we shall assume that the interaction term Q 
takes the form 

Q = SH Pm (8) 

where 5 is a dimensionless quantity which for simplicity is taken to be a 
constant. As a matter of fact the present model with the given form of the 
interaction term was very recently compared to observational data [13J. In 
that work it was found that the allowed range for 5, w, flx,o is the following 

- 0.08 < 5 < 0.03 (9) 

- 1.16 < w < -0.91 (10) 
0.69 < tt x>0 < 0.77 (11) 
while the best-fit parameters were found to be 

5 = -0.03 (12) 
w = -1.02 (13) 
n x ,o = 0.73 (14) 
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where fix,o is the present value of the normalized density of dark energy. The 
normalized densities for matter and dark energy are defined by 



Qm = (15) 

n k2 P x cir\ 
QX = JIP (16) 

and each of them < Qj < 1. Then the first Friedmann equation takes the 
form 

n m + n x = i (17) 

The second Friedmann equation is written in terms of the normalized densi- 
ties 

H = ~H 2 (n m + 1 n x ) (is) 

where 7 = 1 + w. If we define N = lna we can write down the equation of 
motion for Vt m or Q x with respect to N . One finds 

Vl' x = -(1 -n x )(S + 3wn x ) (19) 

and 

n' m = (i-n x )(5 + 3wn x ) (20) 

where a prime denotes differentiation with respect to N. We see that fl' m = 
—Q' x , as it should since Q x + Q m = 1. Since only one of the normalized 
densities is independent, we shall consider the equation of motion for Q x 
and for simplicity we shall drop the index X below. It is obvious that there 
are two critical points, namely 

= 1 (21) 

n., - -I (22) 

To determine the stability of the critical points we have to linearize the 
system. For the first critical point we set Q = 1 + Sfl and we obtain 

5Q' = (5 + 3w)5n (23) 

which means that the critical point is stable for 5 + 3w < and unstable 
for 5 + 3w > 0. In fact according to the values obtained observationally the 
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critical point = 1 is stable. For the second critical point Q* >2 = —S/(3w) 
we set Q = —S/(3w) + 5Q and we obtain 

sn' = -3w(i - n^ 2 )sn (24) 

which means that it is unstable. 

Up to now we have determined the number and stability of the critical 
points of the system. Now we shall focus to the stable critical point and 
obtain the values of r, s corresponding to that. Using the equations of motion 
and the definition for the statefinder parameters we obtain 

a = l + w + i(±.-l) (25) 
9 

r = 1 + - wtt x s (26) 

Notice that for the LCDM model, for which 6 = and w = —1, one obtains 
that s = 0, r = 1 at every instant of time and not only at large times. At the 
critical point fix,* — 1 their values are given by 

s = 1+w (27) 
9 

r = i + - w (i + w ) (28) 
2 

Notice that the dependence on the interaction 5 drops out and the values 
of the statefinder parameters at the stable critical point only depend on the 
dark energy state parameter w. For aw~-l, s~0 and r ~ 1, that is the 
point (s, r) is only slightly different from that corresponding to LCDM. The 
(s — r) plane for this model is shown in Fig. 1, in which the different orbits 
correspond to different values of S — —0.08, —0.03, 0.03. For the numerical 
demands we have let iV = lna go up to iV = 10, which is sufficient for our 
purposes. The statefinders r, s as functions of red-shift are shown in Fig. 3 
and Fig. 4. The present values of r, s are as follows 



w = -1.02, 5 = -0.08 


5(0) = - 


-0.030 


r(0) = 1.100 


w = -1.02, 5 = -0.03 


5(0) = - 


-0.024 


r(0) = 1.079 


w = -1.02, 5 = 0.03 


8(0) = ~ 


-0.016 


r(0) = 1.055 


Table 1: The present values of r, 


s for t 


ie four-dimensional interacting dark energy model. 



Finally, in Fig. 5 we show as a comparison where the curves corresponding 
to quintessence and constant w without an interaction are located on the 
(s — r) plane. 
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2.2 Brane model 



We shall now discuss another dark energy model, this time a brane model. 
In particular we wish to study the DGP model jH] taking into account the 
energy exchange between the brane and the bulk [2]c. Let us first review the 
basic formulae that we shall be using. The model is defined by the action 



S = J (fx^f^j (M 3 R ~ A )+J d 4 xV^h {m 2 R - V) (29) 



plus a matter content both in the bulk and on the brane, where R, R are the 
Ricci scalars of the bulk metric qab and the induced metric h^ u respectively. 
The five- dimensional Planck mass is M, the bulk cosmological constant is 
A/2M 3 <0, the brane tension is V, and the induced-gravity crossover scale is 
r c = m 2 /M 3 . Assuming a flat universe and a perfect fluid on the brane with 
state parameter w, the cosmological equations on the brane read 

p + 3(1 + w)Hp = -T (30) 
H 2 = fj, + 2 7 p + (31) 

i + 2Hh,- X + «^^ = *g- (32) 



a - = p-(l + 3 W ) 1P + p X + 6{1 ; 3whp (33) 

a if) 

where T is the term responsible for the energy exchange between the brane 
and the bulk and it is assumed to have the form T = Ap v [T5] with A, v 
constants. The case A > corresponds to outflow while the case A < 
corresponds to influx. The new parameters f3, j, are related to the old 
ones M, m,V,A by 

w (34) 

(35) 

(36) 
(37) 



A = 


2V 12 

m 2 r 2 


n = 


V 2 
6m 2 r 2 
1 


7 = 




12m 2 


P = 


1 
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and adopting the Randall- Sundrum condition, // = — /3vA. The critical point 
analysis has be done in Ref. [2]c and in the case of dust and influx the results 
can be shown in the table below 





v < 3/2 


v = 3/2 


v > 3/2 


No. of F.R 


1 


or 1 


1 


Nature 


A 


A 


S 



Table 2: The fixed points for w=0, influx 

where the second row shows the number of critical points, while the third row 
shows the nature of the fixed points, attractor (A) or saddle (S). Below we 
shall be interested in the v < 3/2 case, since this is when the system possesses 
always a unique attractor. For the critical point analysis it is useful to define 
appropriate dimensionless quantities. Defining 

2 7 p (fy 7 = H 

D 2 ' D 



UJr. 



D 5 



Coty 



(38) 



where D = y^H 2 —/j,, we obtain the equations 



U). 



l+3w)(uj m —T)Z 



m + <^V> — 1 

A 



u-l 



'1-Z< 



Z' = (l-Z' 



-2Z(1 



(i-z'- 



-z 2 
, 1 



1 



-Z 2 -3(l-3w)/3V 1 ^r, 



1 — &m 

Z 2 -3(l-3w)/?VV, 



1— UJr, 



l + 3w 



-UJr. 



with 



■-d/d,T = D 1 d/dt, while the deceleration parameter is given by 

cj™ -Z 2 - 3( 1 - 3w)/3 



q- 



1 


rl+3w 


"z 2 




L 2 



u m —(l—Z )- 



1— UJr, 



(39) 



(40) 



(41) 



(42) 



Now we shall determine the statefinder parameters in terms of Z, u m using 
their definition and the cosmological equations. The expressions are lengthy, 
therefore we choose to express them in term of the deceleration parameter 



q + 2q 2 
r — 1 



z 



(43) 
(44) 
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One understands that by taking the derivative of q with respect to r and using 
the equations of motion for Z, u m , it is possible to express r, s in terms of 
the latter variables. In practice, we solve the system numerically and obtain 
Z, uj m as functions of r. Then we compute q, q' and finally the statefinder 
parameters cj m (r), Z(r). At the critical points, q' — and in this particular 
model it turns out [2]c that q* — — 1 always holds. Therefore we find that 



at the critical points. These are the values corresponding to the LCDM 
model. The (s — r) plane for this model is shown in Fig. 2, in which the 
different orbits correspond to different values ofi/ = l,l.l,1.2. For our com- 
putational needs we have allowed the mathematical time r go up to r = 100, 
although a few tens would also be sufficient. The basic idea pursued in [2]c 
was that our universe is close to the unique attractor of the system studied. 
Therefore today's values of r, s are approximately the ones corresponding to 
the fixed point. Contrary to the previous four- dimensional interacting model, 
we cannot show for a comparison in the same (s — r) plot the curves corre- 
sponding to the brane model and to quintessence and constant w without an 
interaction. The reason is because the curves for quintessence and constant 
w are located close to the cosmological constant point (0, 1) and cannot be 
seen in the plot. Finally, in Fig. 6,7 we show the statefinder parameters s,r 
as functions of red-shift z for three different values of v = 1, 1.1, 1.2. 

3 Conclusions 

To summarize our work, we have considered two dark energy models. In 
the first model we have considered in four-dimensional General Relativity 
a two fluid model in which an interaction is allowed between dark energy 
and dark matter. We have treated dark energy as a hydro dynamical fluid 
with a constant state parameter w and we have not relied on a concrete 
particle physics model for dark energy. Using dynamical system methods 
we have determined the number and stability of the critical points of the 
system. We have found that there is a stable and an unstable critical point 
and have computed the values of the pair (s, r) corresponding to the stable 
critical point. It turns out that these values depend on the dark energy state 
parameter w only and that the dependence on the interaction 5 drops out. 



r* = 1 
s* = 



(45) 
(46) 
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Our conclusion is that for the observational value of w ~ —1, the values 
of (s,r) for the system under study is only slightly different from that of 
LCDM. In the other model we have considered the DGP brane model taking 
into account a possible energy exchange between the brane and the bulk. The 
assumed form for this energy exchange as well as the number and nature of 
the critical points for this model have been previously discussed. Here we 
have determined the values of the statefmder parameters corresponding to 
the unique attractor of the model and we have shown that are identical to 
the values corresponding to the LCDM model. Finally, we have generated 
plots in which we show s,r as functions of red-shift, and the (s — r) plane 
for each model. 
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Figure 1: The (s — r) plane for the four- dimensional interacting dark en- 
ergy model. The different orbits correspond to different values of 5 = 
-0.08,-0.03,0.03. The marked point at (0,1) corresponds to the LCDM 
model. 



r 




s 



Figure 2: The (s — r) plane for the DGP brane model. The different orbits 
correspond to different values of v — 1, 1.1, 1.2. The marked point at (0,1) 
corresponds to the LCDM model. 
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Figure 3: Parameter s as a function of red-shift z for the interacting model. 
The solid line corresponds to 5 = —0.08, the dashed line corresponds to 
5 = —0.03 and the dotted line corresponds to 5 = 0.03. 
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Figure 4: Parameter r as a function of red-shift z for the interacting model. 
The solid line corresponds to 5 = —0.08, the dashed line corresponds to 
5 = —0.03 and the dotted line corresponds to 5 = 0.03. 
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0.05 0.1 0.15 0.2 0.25 0.3 

Figure 5: The (s — r) plane for the four- dimensional interacting dark energy 
model. Here we also show where the curves corresponding to quintessence 
and constant w without interaction are located. The vertical dashed line 
corresponds to constant w and the dotted line corresponds to quintessence. 
The three curves of Fig. 1 here are shown in the up left corner of the plot, 
close to r = 1. 
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Figure 6: Parameter s as a function of red-shift z for the brane model. The 
solid line corresponds to v — 1, the dashed line corresponds to v — 1.1 and 
the dotted line corresponds to v = 1.2. 
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Figure 7: Parameter r as a function of red-shift z for the brane model. The 
solid line corresponds to v — 1, the dashed line corresponds to v — 1.1 and 
the dotted line corresponds to v — 1.2. 
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